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ABSTRACT 

In this paper we describe a self-contained method for performing the spectral-imaging 
deconvolution of X-ray data on clusters of galaxies observed by the ASCA satellite. Spatially- 
resolved spectral studies of data from this satellite require such a correction because its optics 
redistribute photons over regions which are of comparable size to the angular scales of in- 
terest in clusters. This scattering is a function not only of spatial position but also energy. 
To perform a correction for these effects we employ Maximum-Likelihood deconvolution of 
the image (within energy bands of 1 keV) to determine the spatial redistribution, followed 
by a Monte-Carlo energy reassignment of photon energies with position to determine the 
spectral redistribution. We present tests on simulated cluster data, convolved with the various 
instrumental characteristics and the X-ray background, which show that our methodology can 
successfully recover a variety of intrinsic temperature profiles in typical observational circum- 
stances. In Paper-II we apply our spectral-imaging deconvolution procedure to a large sample 
of galaxy clusters to determine temperature profiles, some of which will be used in subsequent 
mass determinations, presented in Paper-Ill. 

Key words: methods: data analysis - methods: numerical - techniques: image processing 
- X-rays: galaxies - galaxies: intergalactic medium - galaxies: fundamental parameters - 
galaxies: cooling flows 



1 INTRODUCTION 

One of the most important goals in the X-ray study of clusters 
of galaxies is the determination of their mass properties. The ra- 
dial distribution of the total gravitating mass in a cluster can be 
constrained by accurate spatially-resolved observations of the gas 
temperature and density, combined with the simple assumption of 
hydrostatic equilibrium. While the accurate determination of den- 
sity profiles is well within the capabilities of previous instrumen- 
tation (e.g. ROSAT), accurate radial temperature profile constraints 
requires good spectral and spatial resolution. 

The ASCA satellite has the spectral resolution required, but the 
spatial resolution provided by its nested foil mirrors is inadequate 
for the analysis of extended sources, such as galaxy clusters. The 
on-axis half-power diameter of the point-spread function (PSF), 
of approximately Sarcmin, increases and becomes increasingly 
asymmetric with off-axis position. However, the essential problem 
with the PSF is that it varies with energy, such that scattering is 
more severe for harder photons. Consequently, if no correction is 
applied to an ASCA dataset then the determination of the spatially- 
resolved temperature characteristics of a cluster will be incorrect. 
For a truly isothermal cluster, the PSF will produce a temperature 



profil e which appears to rise with increasing radius from the cluster 
core (Takahashi et al. 1995). 



Most of the results presented currently in the literature, that 
attempt to correct for the spatial and en ergy-dependent nature o f 
the ASCA PSF, use a method created by Vlarkevitch et al. (1996 ). 
Their procedure relies on having an accurate prescription for the 
emissivity profile of the cluster - usually a /3-model fit obtained 
from a ROSAT surface-brightness profile. Simulated cluster pho- 
tons are generated according to an initial guess for the cluster spec- 
trum, and spatially distributed according to the emissivity profile. 
The photons are then convolved with the spatially- variable PSF in 
different energy bands, and the model cluster compared with the 
observational data. The source spectrum is then varied at differ- 
ent positions until the model and observational data are statistically 
consistent. 

The advantage of their method is that the spatially-variable 
nature of the PSF is easily incorporated, while the disadvantages 
are that the results may be compromised by the applicability of the 
spectral model and the accuracy of the emissivity profile. The lat- 
ter issue may be important because the ROSAT energy band (0.2 — 
2 keV) is si gnificantly different from that of ASCA [0.5 - 10 keV 
- although [Markevitch et al. (19981) use 2 - lOkeV]. Therefore, 



there may be significant mismatch between the emissivity profile 
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described by ROSAT and that required for the ASCA data analy- 
sis, if the cluster exhibits spatial variations in temperature or other 
properties which affect the spectrum, i.e. a cooling flow cluster will 
have many temperature components and a probable variation in the 
column density across the cluster (although this has to be large to 
affect the spectrum above 2 kcV). 



(Molendi et al. 1999) have presented Beppo-SAX results which 



Using the abov e method, Markevitch et al. (1998) (hereafter 
MFSV) and othe rs QMarkevitch et al. (1996|); S arazin, Wise, & 
Markevitch (1998); ponnelly et al. (1998|); Henric ksen & Marke- 
vitch (1996); Markevitch, Sarazin, & Irwin (1996|); M arkevitch & 
Vikhlinin (1997); [Markevitch (1996[ )] have presented temperature 
profiles for many clusters of galaxies. In a particular MFSV anal- 
ysed a sample of 30 objects and found that most of these clusters 
have temperature profiles that decline with radius, which they pa- 
rameterised with an average polytropic index of 7 = 1.24^q'^2- 

This index is close to, and consistent with at the 2-a bound- 
ary of their uncertainties, the convective instability limit of 7 > 
5/3 = 1.67. While convective instability may be expected in clus- 
ters which have been disturbed, this should not be the case for re- 
laxed clusters. Indeed, because cooling flows should be disrupted 
in significant mer ger events they can be t aken as an i ndicator of 
cluster relaxation (Buote & Tsai 1996 and Buote 1998). However, 



as MFSV themselves find that approximately 60 percent of their 
clusters contain a cooling flow, there is a inconsistency between 
the apparent proportion of relaxed clusters in their sample and the 
possibility that the sample is close to convective instability - as 
implied by their steep temperature gradients. 

Additionally, the steep t emperature gradien ts exacerbate the 
baryon proble m in clusters (|White et al. I993|; Briel, Henry, & 
Bohringer 1992; [White & Fabian 1995[ ) - i.e. the apparent dis- 
crepancy between the large relative fraction of the total mass in 
baryons (i.e. essentially that seen as X-ray emitting gas), and the 
fraction expected from the production of baryons during primor- 
dial nucleosynthesis in a flat (Qo = 1) Universe. Declining tem- 
perature gradients imply smaller total cluster masses than would 
be inferred from isothermal temperature profiles, (if the electro n 
and ion temperatures are in equipartition, Ettori & Fabian 1998 ), 
thereby leading to even larger baryon fraction determinations. 



Although Markevitch et al themselves show that their aver- 
age temperatures determinations agree well with those obtained 
from previous broad-beam satellites such as Einstein, EXOSAT and 
Ginga, the crucial question is whether the shape of their tempera- 
ture profiles are correct. For this their supporting evidence relies 
mainly on the consisten cy of the temperature profiles fr om the dif- 
ferent ASCA detectors. {Markevitch & Vikhlinin (1997 ) have per- 
formed a comparison of ASCA results with ROSAT temperature 
profiles, and found reasonable agreement, but the ROSAT tempera- 
ture determinations have rather large uncertainties. 

There are a growing number of results p ublished which use 
differe nt methods. Many of these results (i.e . [ikebe et al. 1997 



A780; ^ujita et al. I996|: A 399 and A40I; ^zawa et al. 1997 



AWM7- phashi et al. 199"7| : 3A0336+098, MKW3s, A1795 and 
PKS2354— 35) have been compared with those from Markevitch 



et al. by Irwin, Bregman, & Evrard (1998). They highlighted the 



fact that many of these other authors determine that these clusters 
have isothermal temperature p rofiles, even in cases where Marke- 

also presented their own 



vitch et al. find a clear decline. Irwin et al. 



show that A23I9 is isothermal, although these authors are reluc- 
tant to claim any significant discrepancy with the Markevitch et al. 
results on the basis of this one observation. 

Given the importance of generic temperature declines in clus- 
ters, and their contradiction with other results, it is essential that 
the ASCA data are analysed by independent means to check the 
MFSV results. Also, ASCA data currently still provide the best op- 
portunity for the most accurate temperature profiles determinations 
for a large sample of clusters. Thus, we we have created our own 
'spectral-imaging' deconvolution (SID) procedure. Our method is 
self-contained: it requires only ASCA observational data on the 
cluster, background, and the energy-dependent PSF. It does not 
require ROSAT constraints on the emissivity profile. Also, no as- 
sumption about the various components contributing to the source 
spectrum is required in order to perform the deconvolution, and our 
method is essentially non-parametric. The main assumption that we 
have made is that a spatially invariant PSF is sufficient for our pur- 
poses. This has been done for computational speed and for conve- 
nience (we use a third party image deconvolution routine - see Sec- 
tion 2.1), although the consequences of this assumption are tested 
for, and are shown to be acceptable. 

Such tests are central to our methodology, in order that we may 
be confident in the validity of our results. Thus, we have tested our 
procedure on simulated clusters by subjecting them to the various 
instrumental and background effects that will occur in the practical 
application of the method to real GIS data. For example, we have 
investigated the effect of the position-dependent nature of the PSF, 
looked at the degradation in performance with decreasing signal- 
to-noise, and tested the procedures ability to recover a variety of 
intrinsic cluster temperature profiles. 

In summary, this paper presents a spectral-imaging deconvolu- 
tion method for practical use w ith ASCA data on clusters of galax- 
ies. In Paper-II ( White 1999a ), we apply this method to a large 
sample of ASCA GIS cluster observations, and compare our results 
with the results from MFSV. 



2 METHOD 

Our procedure for the spectral-image deconvolution (SID) of ASCA 
data can be broadly divided into two sections: image deconvolu- 
tion followed by spectral reassignment, as summarised by a flow- 
diagram of the method shown in Fig.|l|. 



2.1 Image Deconvolution 

The image deconvolution is performed in the Interactive Data Lan- 
guage (IDL) environment using a maximum-likelihood (M-L) pro- 
cedure, distributed with the ASTROnQ library of contributed rou- 
tines by W. Landsman. Although this procedure accounts for the 
effect of Poissionian noise in the image, it has two significant draw- 
backs: it assumes, (i) that the PSF is spatially invariant, and (ii) that 
the data are monochromatic. 

The first issue of the spatial variance of the PSF is neglected. 



analysis of ROSAT colour profiles which showed that their clus- 
ters were generally consistent with isothermality. Most recently. 
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as we will assume that the PSF appropriate to the position of the 
cluster in the detector image can be used for the whole image de- 
convolution. We argue that this is reasonable, on the basis that most 
of the photons we are interested in arise from the core of the clus- 
ter. However, the i mpact of this assumption is assessed in the test 
detailed in Section 4.5.2. 



The monochromatic issue is dealt with by dividing the data 
according energy, and executing the deconvolution on the image 
for each energy band. This can not be done for the data in each 
individual PI (pulse invariant) channel, as the M-L deconvolution 
procedure requires more counts to work with than will be found 
within one PI channel, and the PSFs which are available in the 
ASCA calibration database (CALDB) at HEASARC (High-Energy 
Astrophysics Archive, Goddard Space Flight Center) are only sup- 
plied in 1 keV bandwidths from 1 to 10 keV. 

Our methodology is to take the X-ray events from each cluster 
observation and divide them into broad spectral energy bands, cor- 
responding to the energy ranges of the PSF images. We then con- 
struct images from the data events in each of these energy bands 
(Figjl] - steps D.l and D.2) and then pass them, along with the 
appropriate point-spread function, to the image deconvolution pro- 
cedure (step D.3). After 50 iterations the M-L procedure returns the 
deconvolved image for that energy band (step D.4). (The effect of 
varying the number of iterations is investigated in Section 4.5.) 



These deconvolved images are then supplied to the spectral 
resassignment phase to determine the PI energy of events which 
will comprise the deconvolved dataset. 



2.2 Spectral Re-assignment 



To make the correct reassignment we need to choose a 'prior' 
spatial distribution for the events within each 1 keV energy band. 
For this we use the M-L deconvolved image (while MFSV use the 
ROSAT data to obtain a higher spatial-resolution emissivity pro- 
file). We can then effectively 'ray-trace' events from this image 
through the telescope optics and determine the most likely asso- 
ciation between events in the deconvolved and convolved planes. 
As we know the energies of each convolved event we then know 
the energy of the deconvolved event. 

In our ray-tracing method we assume the PSF is invariant 
within each energy band. (We discuss the systematic bias that this 
can introduce into our results in Section 4.5.) For each PI chan- 



nel we randomly generate the appropriate number of jowrce-plane 
events from the deconvolved image (steps R. 1 to R.3). We then use 
the PSF to scatter these events into the ;>!ja^e-plane (step R.4). If 
any of these /mage-plane events correspond to the positions of an 
observed photon (for the PI energy in question), then we consider 
this 'mapping' to be successful (R.5). We do this until all the pho- 
tons in the PI channel in question have been successfully mapped. 

Having obtained the position and energy of a i'OMrce-plane 
event we can store the information (R.6) and eliminate the appropri- 
ate photon from the observed event list. We continue the ray-tracing 
until all the observed events are accounted for (R.7). As the number 
of events to be reassigned declines, the probability of finding a suc- 
cessful mapping reduces. For efficiency we mask out parts of the 
PSF which can not provide a successful link between a source- and 
image-plane position. Each randomisation is thereby guaranteed a 
successful mapping and the computational efficiency is improved 
significantly. 



From the M-L procedure we know the spatial probability distribu- 
tion of photons in the deconvolved plane. We also know that the 
total number of counts are conserved between the convolved and 
deconvolved planes. Therefore, we can create a random set of po- 
sitions and energies from the deconvolved image and the observed 
PI distribution. However, we do not know which event has which 
PI energy. 

Instead of attempting to recover the individual PI information 
for each event we could have performed a colour analysis of the 
deconvolved surface-brightness profiles in each 1 keV energy band 
(approximately 85 PI channels). However, this would fail to utilise 
the superior spectral capabilities of ASCA, and would lead to a de- 
generacy in temperature and abundance determinations. Thus, we 
have endeavoured to maintain the spectral resolution of ASCA by 
finding a way of assigning energies to the events list created from 
the deconvolved images. 

The obvious procedure would be to randomly assign an en- 
ergy for each event using the overall spectral distribution of events 
within the 1 keV energy band, however this is not correct if there 
is spatial variation of the PI events over this energy band, as the 
following example illustrates. Consider two spatially distinct point 
sources within a single 1 keV energy band interval: one emitting 
1.2 keV photons and the other 1.8 keV events. A random assign- 
ment of energies in this energy band, which does not take into ac- 
count any spatial information, will result in each point source hav- 
ing half of the 1.2 keV photons and half of the 1.8 keV photons. If 
the same number of photons were detected from each source, they 
would both end up with an incorrect mean energy of 1.5 keV. 



2.3 Practical Considerations 

We have implemented our procedure in the Interactive Data Lan- 
guage (IDL) environment. The deconvolved events listings are writ- 
ten to standard (FITS) format files, which are then analysed using 
conventional X-ray data processing procedures (FTOOLS 4.2 and 
XSPEC 10.00 - see below). We currently only analyse GIS data as 
it has a larger field-of-view (FOV)|^ than the SIS - which allows 
us to determine cluster properties to larger radii. (The larger FOV 
also conserves more of the photons scattered to larger off-axis an- 
gles by the PSF. The SIS also has gaps between the CCDs, and the 
individual CCDs have differing instrumental gains.) 

With the methodology we have defined above, there are vari- 
ous assumptions and complications which have to be assessed be- 
fore we can be confident in applying it to real data. Firstly, we 
have noted, the M-L image deconvolution procedure only accepts a 
spatially invariant point-spread function, while the X-ray telescope 
(XRT) of ASCA exhibits a significant variation with position around 
the field of view (e.g. see Fig. ^). In the spectral reassignment we 
also assume that the PSF does not change significantly across each 
1 keV energy band. Secondly, observational data also include con- 
taminating events from cosmic X-ray sources, and from the X-ray 
detector itself. Finally, systematic effects introduced by the proce- 
dure have to be minimised. 



t The usable area of the GIS is circular and ^ 40 arcmin in diameter, 
while the SIS is square and ~ 18 arcmin on a side. 
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Figure 1. Flowchart of the spectral-unagmg deconvolution procedure. 'SRC refers to the 'source plane' data (i.e. the actual data before it is 
affected by the satelUte instrumentation) which we are attempting to recover, and 'IMG' refers to the 'image plane' data (i.e. the observed data 
after convolution etc.). 
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Figure 2. 

In plots (a) and (b) we show (in logarithmic intensity) the combined XRT and GIS PSF (1 — 2kcV) at two extreme positions: l.Sarcmin olf-axis with 
a position-angle of 5° and 13 arcmin off-axis with a position-angle 40°, respectively. This highlights the change of the PSF with position in the detector 
image. 



Dealing with the last of these points first, our experience has 
shown that a single realisation of the deconvolved data is subject to 
systematic effects. In particular, we noticed that the abundance pro- 
files show abrupt variations which were clearly due to systematics 
in the procedure (the M-L image deconvolution). However, these 
can be reduced by repeating the SID procedure many times. Due to 
time constraints, we found that 10 repetitions was the maximum 
number we could deal with practically. To perform the Monte- 
Carlo procedure we randomise each input (convolved) dataset, and 
then process each deconvolved events list independently in separate 
spectral analyses. The spectral results are then combined to give 
the average profiles presented in the simulation results plots (this 
radial-profile averaging procedure is described in Paper-II). As the 
error bars in these plots are larger than the variations from bin-to- 
bin (and larger than the errors on either the original or convolved 
data) it is clear that the errors in the deconvolved profiles are dom- 
inated by systematic errors rather than statistical errors. However, 
the abundance profiles, in particular, are now much better behaved 
than found in any single realisation. 

For the other two points listed above, we believe the only ef- 
fective way to investigate the validity of our procedure, especially 
in circumstances typical of an ASCA observation, is to create simu- 
lated cluster datasets and then apply the various observational and 
systematic effects to see how well the original physical properties 
of the intracluster gas are recovered by the deconvolution proce- 
dure. In the description below, we indicate how we create these 
various simulations and how we assess the impact of these various 
contaminants in practical observational circumstances. 



3 ANALYSIS 



3.1 Generating Cluster Models 

We generate the basic simulated cluster dataset from a radial den- 
sity profile, parameterised by a standard /3-model, and an isother- 
mal temperature profile. The expected flux from each component (a 
second component is added to mimic a cooling flow cluster) of the 
intracluster medium (ICM) is determined by projecting the volume 
emissivity along the line-of-sight and convolving this with the spec- 
tral response of the detector. The distribution of counts in each PI 
detector channel is then obtained by randomly sampling the intrin- 
sic cluster spectrum. This is repeated until the total observed flux is 
obtained for each thermal component. Similarly, the detected spa- 
tial position of each event is determined by random generation from 
the projected emissivity, followed by convolution with the spatial 
PSF. Note, the effect of vignetting is implicitly incorporated by our 
use of a spatially variable PSF to create the simulated datasets, as 
the normalisation of the PSF describes the efficiency of the detector 
at that position. 

The cosmic background (from observations held at the ASCA 
CALDB at HEASARC) can also be included by adding in the ex- 
pected number of background events from an observation of a blank 
area of real sky, and scaling for the relative exposure time of the 
simulated dataset. Because the background data are taken from ac- 
tual observations of blank-sky fields these data are already con- 
volved with the real spatially-variable PSF and include all the var- 
ious sources of noise, such as the GIS calibration source, e.g. see 
Fig. §(b). 

The convolved simulated dataset, incorporating optional in- 
strumental and background effects, is then supplied to our SID pro- 
cedure, along with the PSFs. The PSFs in each 1 keV energy band 
are selected to most closely correspond to the position of the cen- 
ter of the cluster in the detector image. In the following section, we 
compare the results from the various tests which are designed to as- 
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Figure 3. 



This figure presents a compaiison of tlie original (medium-dark, solid line), deconvolved (dark, solid line) and convolved (light, dashed line) surface- 
brightness profiles for the (TEST-NOBGD-3()K) simulations with no background contamination. This enables a direct comparison to be made of how well 
the SID procedure recovers the true surface-brightness profiles in each energyband. Although the SID procedure introduces extra noise, the true surface- 
brightness profiles are clearly recovered well (the deconvolved data are truncated outside 20 arcmin), and are significantly better representations than the 
surface-brightness profiles from the convolved data, even in the higher-energy bands which have few counts. 



sess the impact of the assumptions imphcit in our method, and the 
various observational contaminants such as the X-ray background. 

We judge the success of the SID procedure by its ability to re- 
cover the properties of the original simulated cluster, and to do this 
we perform a spatially-resolved spectral analysis on the original, 
convolved, and deconvolved events and then compare the results 
graphically. Although the correct recovery of the original tempera- 
ture profiles by the SID procedure implicitly requires that the orig- 



inal surface-brightness profiles are correctly recovered by the M-L 
deconvolution in each 1 keV energyband, this is explicitly shown 
in Fig. |. 
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Figure 4. 

This figure compares the spatial and spectral characteristics of the original, 
convolved and deconvolved datasets. In (a) we show the image of the (TEST- 
GCF-30K) events before spatial convolution with the PSF and the addition 
of background events. The inclusion of these effects produces in the image 
shown in (b) - the cahbration source in the top-right of the image comes 
from the blank-sky background data. The image shown in (c) is that obtained 
after the spectral-image deconvolution analysis - it is truncated at 20 arcmin 
radius from the centre of the field (before the deconvolution) to eliminate the 
calibration source in (b). (These images have a pixel scale of approximately 
0.25 arcmin, have been smoothed with a Gaussian of 2-pixel width, and are 
displayed in logarithmic intensity). Plot (d) shows the (1 — 9 keV) azimuthal 
profiles of the deconvolved (solid line), convolved (dashed), and the original 
data (dotted line) (from non-smoothed data). The upturn in the convolved and 
deconvolved profiles is due to the background, which is not included in the 
original data. In the final panel (e), we compare the spectral PI distribution of 
events extracted from within 3 arcmin of the peak position. The lower line 
is the spectrum from the convolved dataset. The distribution of the original 
and deconvolved events, from the same aperture, are plotted above as a black 
line and lighter line respectively. This plot shows, firstly, that the effect of 
the PSF is significant - the spectral distribution of the convolved events from 
the central 3 arcmin are lower in normalisation (emissivity) and also exhibits 
a different slope (temperature). Secondly, the plot shows that the spectral- 
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3.2 Spectral Analysis 

The spectral analysis procedure which we employ is fairly standard. 
We start by defining annular regions of interest, centred on the peak 
of the X-ray emission, and then fit a spectral model to each of these 
regions and compare the results. The maximum radius allowed for 
any cluster extraction is 20 arcmin from the centre of the field of 
view (although the region of most practical interest will generally 
be within the inner 10 arcmin). This eliminates the gain uncertain- 
ties at the edge of the GIS detector and contamination from the 
calibration source. However, if the background has been incorpo- 
rated the maximum radius is usually reduced to the radius where the 
background-subtracted surface-brightness profile remains positive. 
The background spectrum is extracted from the blank sky observa- 
tions using the same spatial regions defined by the cluster data. In 
this way any gross position-dependent detector variations, such as 
instrument gain, should cancel out. 

After defining each annular region, such that each annulus 
contains a certain fraction of the total (if applicable - background 
subtracted) counts (0. l|^ - which gives a theoretical maximum of 
10 radial bins; see Paper-II for details on the methodology), we ex- 
tract events using FTOOLS 4.2. This is done for both the source 
and background spectra. These spectra, together with the ancil- 
lary response matrix {i.e. ARF file; determined using ASCAARF), 
are passed to the XSPEC 10.00 spectral analysis package ( Arnaud 
1996), and fitted with the chosen spectral model. We note that, until 
this point, no assumption has been made about the intrinsic cluster 
spectrum, and any chosen model can be fitted to the deconvolved 
data - contrary to the procedure used by MFSV which requires a 
spectral model to be defined as an integral part of the deconvolu- 
tion. 

In this analysis we fit a single temperature plasma emission 
modej^ - regardless of whether the simulated data were created 
with two thermal components. This means that we are not using the 
correct model for the data, but for the simulations this is not impor- 
tant. The crucial fact is whether we recover the same results from 
the deconvolved data as from the original data using the same spec- 
tral model. Thus, in the following results we are looking for good 
agreement between the spectral fits to the deconvolved and origi- 
nal datasets (i.e. simulated data which have been convolved with 
the spectral response, but not the spatial PSF). (We do note, how- 
ever, that our lack of a cooling flow component, and particularly 
the excess absorption, in the fitted model is largely responsible for 
the radial increase in metallicities in the cooling flow simulations, 
which is not apparent in the single thermal component models.) 

Finally, we note that if background contamination is included 
in the test it is also applied to the convolved data (and thereby the 
deconvolved data), but not the original data. Thus, no background 
subtraction is performed in the spectral analysis of the original 
dataset. This also means that (except in the test which has no back- 



t The nominal fraction of the total number of background-subtracted 
counts per annulus is 0.1, the soft-limit on the minimum number of counts 
in an annulus is 1,000 and the hard-limit is 500 counts. 
^ Our specific model is that of a single MEK AL thermal component {i.e. 
Mewe, Gronenschild. & van den Oord 19851 Mewe. Lemen, & van den 



TEST-TAK-30K 
Rodius (Mpc) 




Rodius (orcmin) 



Figure 5. This figure shows the effect of the PSF on isothermal cluster data 
(the underlying original isothermal data is shown by the square symbols with 
the dotted line). Using a model similar to that of Takahashi et al. (which is 
shown in the temperature profile panel only - by the square symbols with 
dashed-line error bars) we also obtain the radially increasing temperature pro- 
file in the uncorrected convolved data (dashed line, triangle symbols). This last 
profile (solid line with diamond symbols) shows the recovery of the isother- 
mal profile from the convolved data using our SID procedure. (See Fig. ^for 
other plot details.) 



bord 1986), absorbed (|m orrison & McCammon 1983| ) by a foreground col- 



umn density of foreground material. 
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ground added) the convolved and deconvolved results are subject 
to extra noise which is not in the original dataset spectral fits. 



4 CLUSTER SIMULATION RESULTS 

4.1 Takahashi et al. Model - The Effect Of The PSF On 
ASCA Cluster Data 



TEST-GCF-30K 
Rodius (Upc) 



In the Introduction we referred to a paper by ( Takahashi et al. 1995 ) 



which showed that the energy-dependent PSF of ASCA can give 
rise to an apparently increasing temperature profile in an isother- 
mal cluster, if the data were not corrected for the effect of the PSF. 
Before we discuss our various simulation results, where we inves- 
tigate different observational and physical scenarios, we show that 
we can reproduce this systematic effect when using similar param- 
eters. (The magnitude of the effect depends on parameters, such as 
apparent size of the cluster core.) 

As the effect of the PSF becomes severer for hotter clusters, 
we chose to approxima te the hottest temperature profile modeled by 
Takahashi et al. (1995). We do this using a single /3-model with a 
velocity dispersion of 1, lOOkms"^ (which with the /3spcc = 0.8 
which we use in our other models gives a central temperature of 
approximately 9 — lOkeV), an index for the profile of /Sat = 0.6 
and a core radius of -Rcorc = .13 Mpc (resulting in a 1 arcmin 



core, similar to Takahashi et al. ). Further details of the parameters 
used are given in Table |l| under the TEST-GCF-30K model. 

In Fig. ^ the effect of the PSF can be seen in the convolved 
dataset's temperature profile (dashed line with triangle symbols), 
where the temperature in the core is low and the n it gradually in- 
creases with radius. This trend is similar to the Takahashi et al 



profile, which we present in the plot as square symbols with dashed 
error-bar lines. The diagram also shows that our SID method recov- 
ers the true radial temperature, metallicity and emissivity profiles 
from the convolved data. 

In the following tests we create a suite of models to investi- 
gate the effects of various observational and physical conditions on 
the ability of the SID procedure to recover the true intracluster gas 
properties. 



4.2 Standard Test - A Cluster With A Strong Core 
Temperature Decline 




Rodius (orcmin) 



With our ba sic p rescription for creating a simulated cluster dataset 
(see SectionBjj), we choose our standard model to represent a 'Gi- 
ant Cooling Flow' (TEST-GCF-30K). The variation in the average 
temperature in this type of cluster is expected to be significant, and 
will therefore provide a severe test for the SID procedure. To pro- 
duce a core temperature decline we add a second cooler thermal 
component to the core region. Although this is not a physically 
realistic representation of a cooling flow it mimics a core temper- 
ature drop (albeit rather sharpl y) which is seen in ty pical cooling 
flow clusters (e.g. Centaurus - Allen & Fabian 1994 ). (Remember, 



we are not necessarily interested in the physical correctness of our 
model but whether the SID recovers the intrinsic energy and spatial 
distribution of the simulated cluster - whatever that might be.) 



Figure 6. In this figure we plot the radial variations of electron tempera- 
ture, metallicity (fraction of Solar), and model noiTnalisation {i.e. emission 
measure: 10^''/(47rD|^ ) J ricn-iidV) for a single thermal (MEKAL) com- 
ponent fit to the standard model {TEST-GCF-30K; see Section^). The results 
for the deconvolved data are shown by the soHd line, the convolved data by 
the dashed line, and the original data by the dotted line. The error bars on the 
temperature, abundance and nonnalisation are one standard deviation, while 
those on the x-axis data show the radial extent. 



The parameters have been chosen to approximate the physical 
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Test Name 


Comments 


Systematic/Instrumental Flags 




Cluster Component 




Fig. 
Ref. 


Bgd. 


Pos. 


Exp. 






ATVh " 




R " 

^corc 




TEST-TAK-30K 


Takahashi et al. model 


Yes 


4.5 


30 


50 


9.87 


N/A 


N/A 


N/A 


N/A 1 


q 


TEST-GCF-30K 


'Standard' model 


Yes 


4.5 


30 


50 


6.61 


6.8 


1.65 


0.2 


r 

0.3 1 




TEST-MCF-30K 


Smaller cooling flow effect 














3.31 






|a) 


TEST-ISO-30K 


Isothermal 












N/A 


N/A 


N/A 


N/A 


5(b) 


TEST-DTX-30K 


Temperature decline 










1.65 





6.61 


0.8 


5 


|c) 


TEST-NOBGD-30K 


No background 


No 




















TEST-GCF-15K 


Short exposure 






15 














(a) 


TEST-GCF-60K 


Long exposure 






60 














If) 


TEST-NIT25-30K 


Less image deconv. iterations 








25 














TEST-NIT150-30K 


More image deconv. iterations 








150 












m) 


TEST-OFF-30K 


Off-axis cluster 




9.0 
















J(b) 



Table 1. Variable Parameters For Generating Simulated Clusters 

In all but the single-phase temperature model (TEST-ISO-30), there are two gas phases, as indicated by the primary and secondary superscripts. 
The model components which are inapplicable (i.e. the TEST-ISO-30K test) are indicated as such by 'N/A'. Entries with '..' indicate the 
values default back to the 'standard model': TEST-GCF-30K. Parameters which are invariant between each test are: (i) the total column 
density (i.e. Galactic; acting on both components where applicable): 1.36 X 10^^ cm^^ , (ii) redshift: 0.0881, (iii) baryon fraction: 0.1, (iv) /3- 
parameter: 0.67, (v) abundance: O.4Z0, (vi) velocity dispersion: 900 kms^^ (both components where apphcable), (vii) core radius: 0.2 Mpc 
(primary component), and (viii) maximum radial extent: 5 Mpc (primary component). Each test is labelled to indicate the primary point in 
question. Those tests which investigate different cluster temperature profiles are: giant cooling flow: 'TEST-GCF-30K' ; moderate cooling 
flow: 'TEST-MCF-30K' ; isothermal cluster: 'TEST-ISO-30K' ; radially decreasing temperature: 'TEST-DTX-30K\ Those tests investigating 
instrumental/systematic effects are: inclusion of cosmic background 'TEST-NOBGD-30K' ; off-axis position 'TEST-OFF-30K'; number of 
image-deconvolution iterations: 'TEST-NIT25-30K' and 'TEST-NIT150-30K'; exposure time (in kilo-seconds): 'TEST-GCF-15K' and 'TEST- 
GCF-60K' . The abbreviations for the the systematic/instrument flags are as follows: 'Bgd.': background included; 'Pos.': off-axis position of 
the centre of the simulated cluster; 'Exp.': exposure duration; and 'Njt': number of iterations used in the M-L image deconvolution stage. 
The abbreviations for the cluster properties are: gas temperature 'Tx'; excess absoiption (only on second component): 'AA'^h core-radius: 
'RcoTc -land truncation maximum radius: '/Jmax'. Finally 'Fig.' indicates the figure reference for each test, while each test is discussed in 
Section W. 



characteristics of Abell 478. This is a cluster of moderate temper- 
ature, and distance, but is fairly luminous due to its large cooling 
flow. It should be remembered that some of the conclusions drawn 
from the following tests will be dependent on the characteristics 
of the cluster in question, e.g. whether the cluster has a cooling 
flow, or not, and how bright the cluster is, etcetera. Table |l| shows 
the parameters we have use in each test. The primary component 
has a temperature of 6.6 keV; the cooler component is at 1.65 keV, 
but its radial extent is truncated at 1.5 x i?corc = 0.3 Mpc to 
limit it to the centre of the cluster. The cooler component spec- 
trum is also modified by an additional column of absorbing ma- 
terial of five times the assumed Galactic value (which is set at 
1.36 X 10^^ cm~^ ) to model excess absorption (i.e. W hite et al. 
1991 and [lohnstone et al. 1992| ). For this standard model, the ob- 
servational characteristics are an exposure of 30 k s, with the cen- 
tre of the cluster placed 4.5 arcmin off axis (fairly standard for an 
ASCA GIS observation), and the data are convolved using a a spa- 
tially varying PSF. (The PSF is interpolated from the grid of 11 
different positions for the PSF in each energy band, supplied in the 
CALDB.) 

The effect of the spectral-image deconvolution on the con- 
volved data becomes apparent when we look at the radial varia- 
tions in the spectral fit parameters in Fig. |^ This figure shows that 
the variations of the original ICM properties are recovered by the 
SID, and are a much better representation of the original variations 
than those obtained from the convolved data. (If a hydrostatic mass 
determination were made from the convolved data on such a clus- 
ter, then careful consideration of the complications would be re- 



quired, if meaningful results were to be obtained.) We note, how- 
ever, that the abruptness of the temperature profile decline in the 
original data, which is rather artificial due to the construct of our 
model, is somewhat smoothed over in the deconvolved temperature 
profile. 



4.3 A Variety Of Different Temperature Profiles 

In the next three tests we have attempted to span a wide range of 
possible temperature profiles. With the 'Medium Cooling Flow' 
model, TEST-MCF-30K, we investigate whether a less dramatic 
temperature decline can be recovered. In this case, the simulated 
data were created by adding a secondary component which is 
slightly hotter, i.e. 3.3 keV, than in the previous model. Figure ^(a) 
shows that the small differences between the original and con- 
volved temperature profiles are still resolved with the deconvolved 
data. 

In the TEST-ISO-30K model. Fig. |(b), we eliminate the sec- 
ond cooler component altogether, and produce an isothermal tem- 
perature profile. This test shows that the deconvolution procedure 
essentially has no systematic effects which severely distort the re- 
covery of a flat temperature profile and, for example, results in a 
temperature decline. 

In the last temperature profile variation, TEST-DTX-30K - 
shown in Fig. |7|(c), we have created a temperature profile which 
decreases radially outwards, to approximate the general trend for 



Spectral-Imaging Deconvolution Method 1 1 




(a) (b) 

Figure 7. This figure summarises tlie tests of differing temperature profiles, as follows: (a) moderate core teipperature drop, (b) and isothermal temperature 
profile, and (c) and temperature decline with radius. The line styles etc. are the same as those defined in Fig. pi 



decreasing temperature profiles found by MFSV. This model was 
created by making the hotter component more extended {i.e. the 
core-radius is larger), than the cool component (note, the truncation 
is set at 5 Mpc for both components and there is no excess absorp- 
tion on the cooler emission component). Although the difference 
between the original and convolved temperature determinations in 
any particular radial bin is small, the general slope of the original 
profile is more accurately reproduced by the deconvolved data. 



Despite the lower count rates from the TEST-MCF-30K, 
TEST-ISO-30K models, as they incorporate less flux from the sec- 
ondary temperature component, these results show that the success 
of the deconvolution does not appear to be dependent on the intrin- 
sic temperature variations, and that there are no severe systematic 
biases in the deconvolved profiles. 
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(c) 

Figure 8. This figure summarises the test where no baclcground is included, 
Figure 7. - continued. aid therefore shows that performance of the SID method in optimal circum- 



stances. 



4.4 Statistical Effects Due To Background And Exposure 

We now turn our investigation to the impact of the various instru- 
mental and systematic effects which occur in ASCA observations 
of clusters: the effect of background event contamination, and the 
impact of signal-to-noise variations. 

In TEST-NOBGD-30K we see the intrinsic ability of the SID 
procedure in the absence of any background contamination. The 



results are clearly more accurate, but are still qualitatively similar 
to the results from the standard simulation data which includes the 
background. This shows that performance of the SID procedure is 
relatively unaffected by background noise, at the expected level. 

In Figs. ^ we show the effect of variations in data quality, using 
exposure times of 15 and 60 ks, i.e. factors of a half and twice the 
standard exposure (TEST-GCF-15K and TEST-GCF-60K). These 
models include the cosmic background, and therefore the reduction 
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(a) (b) 

Figure 9. This figure summarises tlie test wliicli investigates tlie effect of signal-to-noise through different exposures [(a): 15 and (b): 60 ks, i.e. a factor of 
one-half and twice the standard test's exposure]. 



in exposure time also results in a degradation of the signal-to-noise 
ratio. 

In both cases the deconvolved temperature and emissivity pro- 
files form much better representations of the original data than the 
convolved data. (The performance will, of course, differ for clus- 
ters with differing count distributions i.e. surface-brightness pro- 
files, and count rates, i.e. signal-to-noise for a given exposure.) 



4.5 Systematic Effects From Our Procedure 



In Section 2.3 we discussed the Monte-Carlo technique which we 



have used to average-out systematic variations which occur in a 
single run of the SID procedure. We now discuss some potential 
sources for these systematic effects. 

The primary assumptions in our procedure are that we can use 
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(a) (b) 

Figure 10. This figure summaiises the tests which utihze fewer (25) or greater (150) number of iterations in the M-L stage of the image deconvolution, 
compared to the standard number of 50. 



a spatially invariant PSF for the image deconvolution phase, and 
that the PSF is invariant across each 1 keV energy band for the 
spectral-reassignment. Before we investigate these issues we will 
look at the effect of varying the number of iterations employed in 
the image deconvolution phase of the SID procedure. 



4.5.1 The Number of Image Deconvolution Iterations 

In the TEST-NIT-25 and TEST-NIT-I50 tests we vary the number 
of iterations employed in the (M-L) image deconvolution stage of 
the SID procedure from the standard number of 50, to 25 and 150 
iterations, respectively. Figures Qa) and [l|(b) essentially show 
that the results are similar regardless of the number of iterations. A 
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greater number of iterations will produce sharper features, and thus 
follow the sharp temperature decline in the cooling flow models, 
but this will also exaggerate spurious features (although this will 
be minimised by the 10 Monte-Carlo runs of the SID procedure). 
For this particular cluster model, we find that 50 iterations appears 
to be optimal for recovering of the core temperature decline without 
introducing severe systematic effects into the abundance profile. 



4.5.2 The Assumption Of A Spatially Invariant PSF 



TEST-0FF-30K 



In Section 2.1 we described the use of the Maximum-Likelihood 



procedure, which requires the use of a spatially invariant PSF. How- 
ever, the real ASCA PSF varies significantly with position (e.g. see 
Fig. ^). Thus, we have created all our simulated data using a spa- 
tially varying PSF (interpolating from the CALDB sample of 11 
different PSF positions), and then deconvolved using a. fixed PSF 
which is appropriate for the position of the centre of the cluster in 
the detector image. All the results presented so far show that this is 
a reasonable assumption which does not severely compromise the 
SID results. This is because most of the cluster photons arise from 
the core region, and therefore the fixed PSF approximation we have 
chosen is acceptable. 

Our assumption will be poorer for more extended clusters, and 
for clusters centred further off-axis because the PSF changes more 
rapidly with increasing displacement from the centre of the detec- 
tor (mainly in symmetry - see Fig. In the TEST-OFF-30K test 
we see the results for a cluster placed near the edge of the field 
of view (9 arcmin from the on-axis position). Although the (abun- 
dance profile in particular) results are worse than for the near on- 
axis results (TEST-GCF-30K) the essential trends of the true emis- 
sivity and temperature profiles are recovered. We emphasise that 
most real cluster observations occur at the on-axis position we have 
selected for our standard TEST-GCF-30K test. 



4.5.3 The Assumption Of An Spatially Invariant PSF Within 
Each 1 Kev Energyband 

The scatter of the PSF continually degrades with increasing pho- 
ton energy. This means that for events at the upper end of the en- 
ergy band we will be using a PSF which is too narrow, and for the 
lower energy data we will be using a PSF which is too wide. When 
the spectral-reassignment procedure determines a map of events 
between the source- and image-plane, it will choose hard events 
which are too concentrated, and soft events which are too far away 
than would be correct. The precise effect of this is complicated and 
probably depends on the spatial distribution of events, but undoubt- 
edly it will change the slope of the spectral distribution within each 
1 keV energy band, resulting in a slight discontinuity between the 
boundaries of each energy band. Fortunately, statistical noise in the 
data will help obscure these discontinuities, and none are visible in 
the deconvolved spectrum of the central region of our test shown in 
Fig. §(e). 

In summary of all these tests, we believe that our procedure 
will provide accurate results under a wide-range of observational 
circumstances. 




Rodius (orcmin) 



Figure 11. This figure shows the limitation of using the fixed-spatial PSF in 
the SID procedure when the cluster is positioned at a large (9 arcmin) off-axis 
angle. 



5 CONCLUSIONS 

In this paper we have presented a method for performing a spectral- 
imaging deconvolution analysis of X-ray data which has been af- 
fected by an energy- and position-dependent point-spread function 
(PSF), such as that produced by the ASCA satellite mirror arrange- 
ment. Without correction any spatially-resolved spectral analysis of 
ASCA data will be incorrect, as higher energy photons are scattered 
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further than lower energy photons. In studying spatially extended 
sources, such as clusters of galaxies, it is essential that a correction 
is made for this effect if accurate conclusion are to be drawn from 
the radial variations in the spectral properties of the object, e.g. the 
temperature and thereby mass profiles of galaxy clusters. 



Markevitch et al. (1996) has presented a procedure which at- 
tem pts to correct for the ASC A PSF, and in an analysis of 30 clus- 
ters Markevitch et al. (1998) found that virtually all their temper- 



ature profile results decrease with radius, according to an average 
polytropic index of 7 = 1.241q'^2- If correct, this result is very 
important as it implies that clusters have less mass than would be 
assumed if they were isothermal, and thereby has implications for 
cosmology through cluster mass distributions and baryon fractions. 
However, the fact that the polytropic index is so steep, and close 
to the convective instability value of 7 = 5/3 — 1.67, is of some 
concern as many clusters in their sample are cooling flows, and 
should therefore be relaxed clusters. It is therefore imperative that 
the MFSV results are verified by independent means. 

We have developed an independent method which is self- 
contained, in that it does not require observational data on the same 
object from other X-ray satellites, or any assumption about the 
source spectrum. This means that potential systematic problems, 
arising from cross-calibration uncertainties, are eliminated. Also, 
unlike Markevitch et al., we present tests of our deconvolution pro- 
cedure on realistic simulations of ASCA observations of clusters to 
show that our method successfully recovers the true radial proper- 
ties of the intracluster gas from data under a wide variety of condi- 
tions. In particular, we have tested the recovery of differing intrinsic 
temperature profiles and we have included contamination from the 
observed cosmic and instrumental background of ASCA, to show 
that the procedure will work in practical application on reasonably 
bright clusters. We have also shown that the various limitations 
(i.e. were are forced to assume a spatially-invariant point-spread 
function) do not unduly compromise the ability of the procedure to 
successfully recover the original temperature profiles. All the de- 
convolved profiles we obtain are in good agreement with the true 
intrinsic properties. 

In conclusion, our spectral-imaging deconvolution procedure 
is successful at recovering the original radial properties of the in- 
tracluster gas - for a range of different physical scenarios, and 
at the expected level of background and instrumental contamina- 
tion. This procedure is, therefore, entirely capable of working in 
practical^ircumstances^i^^ 

1999a) we apply this spectral-imaging deconvolution procedure to 
real A5CA data on 106 clusters, and compare the r esults with those 



from Markevitch et al. 



In Paper-Ill (White 1999b) we select a sub 



sample of these objects for mass profile determinations. 
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